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FOREWORD 


This  report  was  prepared  in  the  Components  Branch  (APTC),  Turbine  Engine  Division,  Air 
Force  Aero  Propulsion  Laboratory,  Wright-Patterson  Air  Force  Base,  Ohio,  under  Project 
3066,  “Gas  Turbine  Technology,”  Task  306603,  “Advanced  Engine  Studies.” with  Charles  E. 
Bentz  as  Project  Engineer. 

This  report  covers  work  conducted  within  the  Components  Branch  in  the  time  period  between 
July  1965  and  June  1967  and  was  submitted  by  the  author  31  August  1967. 


This  technical  report  has  been  reviewed  srd  is  approved. 
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ABSTRACT 


This  report  describes  a  digital  computer  program  titled  SMOTE  (Simulation  of  Turbo  fan 
Engine).  SMOTE  is  a  computer  program  for  balancing-cycle  turbofan  engines  capable  of  run¬ 
ning  both  design  and  off-design  points.  Component  performance  map**  are  reduced  to  Block 
Bata  (tabular  form)  to  preside  a  base  for  nainniating  component  performance.  The  design 
point  i3  run  first  and  map  correction  factor  ^  are  calculated  to  scale  the  components  to  the 
desired  performance.  These  correction  factors  are  then  applied  to  the  component  perfor¬ 
mance  maps  at  off-design  points.  Initially  *  cn  the  program  is  running  at  an  off-designpoint, 
the  cycle  is  not  balanced,  and  errors  (for  ev*mple,  work  required  by  the  compressor  minus 
stork  supplied  by  the  turbine)  are  generated.  Small  changes  in  engine  independent  varianlec 
(for  example,  compressor  speed)  then  prodace  small  changes  in  the  errors,  and  these  dif¬ 
ferential  Changes  are  loaded  into  e  matrix.  The  matrix  is  then  solved  for  the  set  of  independent 
variables  which  results  in  zero  errors,  thus  balancing  the  cycle.  Actually,  this  process  may  be 
repeated  several  times  before  it  reaches  a  balanced  point  because  there  is  a  nonlinear  re¬ 
lationship  between  the  independent  variables  >r=i  the  errors.  Sample  results  are  included  in 
tills  report. 

(Distribution  of  tills  abstract  is  unlimited.) 
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T2 
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bleed  from  fan  lost  to  cycle  (leakage) 

bleed  from  compressor  to  duct  (leakage) 

bleed  from  compressor  to  high  pressure  turbine  (cooling) 

bleed  from  compressor  to  low  pressure  turbine  (cooling) 

bleed  from  compressor  overboard  for  customer  use 

corrected  speed 

delta-H  corrected  for  temperature 
enthalpy 

percent  speed  of  the  compressor 
percent  speed  of  the  fan 
pressure 
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turbine  flow  function,  low  pressure  turbine 
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main  combustor  burning  temperature 
duct-burner  burning  temperature 
afterburner  burning  temperature 
afterburner  fuel  flow 
main-combustor  fuel  flow 
duct -burner  fuel  flow 
gas  flow  rate 

pressure-ratio  ratio  of  the  compressor 
pressure-ratio  ratio  of  the  fan 
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SECTION  I 
INTRODUCTION 

Recent  advances  in  turbine-engine  state  of  the  arthave  increased  the  requirements  for  more 
and  better  cycle  studies.  These  cycle  studies  are  needed  to  monitor  present  engines,  determine 
sensitive  or  critical  areas  in  near  future  engines  now  under  development,  and  to  explore  the 
advantages  and  disadvantages  of  proposed  advanced  engine  cycles  for  future  aircraft. 

Parametric  cycle  Studies,  which  involve  essentially  numerous  design-point  calculations, 
partially  fulfill  this  need,  particularly  for  optimizing  a  cycle  for  a  specific  single  design-point 
mission.  However,  with  multimission  aircraft  being  emphasized  increasingly  and  with  the  need 
for  determining  off-design  performance,  the  requirement  for  a  balancing  cycle  computer 
program  (that  is,  one  which  simulates  a  turbine  «mgine  at  both  design  and  off-design  points) 
becomes  definite  and  essential. 

The  purpose  of  this  report  is  to  describe  a  digital  computer  program  for  balancing-cycle 
turbofan  engines.  The  program,  titled  SMOTE  (Simulation  of  Turbofan  Engine),  was  developed 
in  the  Components  Branch,  Turbine  Engine  Division  Air  Force  Aero  Propulsion  Laboratory, 
to  meet  the  requirements  given  in  the  preceding  paragraphs.  In  addition  to  meeting  these 
requirements ,  SMOTE  is  considerably  more  flexible,  requires  less  computer  storage  or  space, 
and  requires  less  computer  operating  time  than  previous  engine  cycle  decks  of  comparable 
sophistication. 

Part  I  of  this  report  describes  the  method  of  engine  calculations  and  the  balancing  technique 
and  gives  some  sample  results.  Part  II  is  intended  as  a  user’s  manual  and  includes  instructions 
for  setting  up  and  running  the  program,  as  well  as  a  program  listing.  The  parts  may  be  used 
independently  of  one  another. 


AFAPL-TR-67-125 

Parti 

SECTION  n 
SUMMARY 


SMOTE  Is  a  computer  program  for  balancing-cycle  turbofan  engines  which  presently  ises 
component  performarce  maps  for  the  fan,  compressor,  combustor,  and  both  torbir  s  to 
provide  the  basic  performance  data,  but  it  can  easily  be  expanded  to  include  add!,  onal 
component  performance  maps  if  available.  The  maps  are  in  Block  Data  form  and  are  scaled 
Internally  to  simulate  a  specific  engine.  Errors  due  to  an  unbalanced  cycle  are  generated  at  off- 
design  points,  and  the  effects  of  small  changes  in  independent  variables  upon  the  errors  are 
determined.  A  matrix  of  differential  error  equations  is  then  so lred  to  determine  the  correct 
values  of  the  independent  variables  which  would  produce  zero  errors.  A  flow  chart  of  the 
program  is  shown  in  Figure  1. 

For  a  more  accurate  simulation  of  a  particular  engine,  performance  maps  for  other  com¬ 
ponents  could  be  added;  for  example,  duct-burner  or  afterburner  maps  may  be  desired.  It 
should  also  be  mentioned  that  other  formats  for  presenting  maps  may  be  used  as  readily  as 
those  presented  in  this  report.  Rather  than  inputting  bleed  air  values  at  each  point  or  using  a 
constant  bleed,  a  bleed  schedule  could  be  used.  In  addition,  if  a  variable-area  nozzle  is  to  be 
simulated,  a  nozzle  area  schedule  could  be  used.  Or  an  engine  control  system  could  be  used 
which  would  set  fuel  flow,  bleeds,  and  nozzle  areas  as  some  function  of  a  power  lever  angle. 

The  complexity  of  an  engine  cycle  can  be  increased  by  increasing  the  size  of  the  matrix 
(increasing  the  number  of  partial  differential  equations).  For  example,  a  basic  triple-spool 
turbofan  cycle  could  be  represented  using  a  matrix  of  nine  equations.  Or  a  T-compressor  fan 
engine  composed  of  a  fan  tip,  fan  hub,  low  pressure  ccr  pressor  (running  at  the  same  speed  as 
the  far),  high  pressure  compressor,  combustor,  two  turbines,  gas  mixer,  and  afterburner 
could  be  represented  using  a  matrix  of  eight  equations. 
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Figure  1.  SMOTE  Computer  Program  Flow  Chart 
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SECTION  m 
HISTORY 


Until  about  six  years  ago.  most  general  cycle  calculations  in  the  Components  Branch  were 
done  by  hand,  although  some  computer  programs  were  available  for  specific  engines.  About 
that  time  a  turbojet,  parametric-cycle  study  program  (SPEEDY)  was  conceived,  and,  from 
this  program,  a  more  general  turbojet  or  turbofanprogram  (CARPET)  with  many  configuration 
options  was  developed.  CAr~*ET  is  still  in  general  use  for  parametric  and  optimization  studies. 

About  three  years  ago,  a  balancing-cycle  turbojet  or  single-spool  computer  program 
(SSPOOL)  was  developed  within  the  Components  Branch.  The  engine  component  calculations 
were  based  essentially  on  those  »n  CARPET,  and  the  arming  technique,  which  depended  upon 
a  quadratic  interpolation  routine  ( AFQU1R) ,  involved  two  nested  haianmng  loops.  The  inner 
loop  was  balanced  using  PCNC  as  an  independent  variable  and  the  work  difference  between 
compressor  and  turbine  an  die  dependent  variable.  The  outer  loop  was  balanced  using  ZC  (see 
Figure  2  for  a  definition  01  ZC)  and  the  pressure  required  by  the  fixed-area  exhaust  nozzle. 
After  the  inner  loop  was  balanced,  die  outer  loop  was  changed  in  an  attempt  to  it. 

Naturally,  changes  in  die  outer  loop  necessitated  rebalancing  the  inner  loop.  This  method, 
although  rather  crude,  worked  well  for  a  turbojet  cycle. 

The  SSPOOL  concept  was  then  extended  to  a  turbofan  or  dual -spool  cycle  which  resulted  in 
a  new  program  called  DSPOOL.  By  logical  extension  this  required  focr  nested  loops  with  four 
independent  variable*  (PCNF,  ZF»  PCNC,  andZQ  and  four  dependent  variables  or  errors  (two 
work  errors  and  two  nozzle  pressure  errors  for  a  separate  flow  cycle;  or  two  work  errors,  a 
mixing  static  pressure  error,  and  a  nozzle  pressure  error  for  a  mbmH  flow  cycle).  Although 
the  method  worked,  computer  time  was  excessive,  and various  techniques  (such  as  changing  the 
order  of  the  independent  variables  or  using  a  varying  tolerance)  were  tried  in  an  attempt  to 
shorten  the  balancing  time.  These  attempts  were  only  partially  successful. 

Other  balancing  techniques  using  various  mathematical  solutions  were  experimented  with, 
xnd  the  present  method  was  finally  developed.  This  method  involves  no  nested  balancing  loops; 
instead,  a  matrix  is  loaded  with  differential  errors  caused  by  small  changes  in  the  inAepewfemt 
variable*.  The  matrix  is  then  solved  for  the  zero  error  condition.  SMOTE  reduced  computer 
tin.«  by  an  average  factor  of  about  4  as  comnared  to  DSPOOL. 
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SECTION  IV 


METHOD  OF  ENGINE  CALCULATIONS 


1.  COMPONENT  MAPS 

The  performance  of  the  major  engine  components  is  based  on  component  maps.  These  maps 
are  usually  obtained  from  analytical  methods  or  rig-testing  and  are  then  converted  into 
Block  Data  subroutines  for  use  by  SMOTE.  The  maps  presently  included  in  SMOTE  are  eery 
general  and  do  net  represent  any  particular  mgtne  or  engine  components. 

The  component  maps  are  scaled  at  the  engine  design  point  by  SMOTE  in  order  to  match 
their  performance  to  a  desired  set  of  performance  figures  which  are  input  as  data.  Scaling 
or  correction  factors  are  calculated  and  then  applied  to  the  maps  at  off-design  points.  The 
scaling  process  is  ’Jneax;  therefore  correction  factors  near  unity  result  in  the  hipest  accuracy 
of  component  simulation.  This  means  that  the  component  maps  used  should  represent  or  be 
similar  to  the  actual  components  in  the  engine  being  simulated.  However,  with  the  loss  of  a 
little  accuracy,  maps  representing  advanced  components  could  be  interchanged  to  determine 
the  effect  on  the  overall  cycle. 

SMOTE  presently  fncludea  component  maps  far  the  frn,  compressor,  combustor,  and 'both 
turbines.  Duct  burning,  duct  losses,  gas  miring,  afterburning,  tailpipe  losses,  and  nozzle 
losses  are  all  calculated  or  input,  but  these  characteristics  could  also  be  included  as  Block 
Data  if  maps  were  available.  Likewise,  schedules  far  bleed  sir  and  variable  area  nonles  could 
be  used. 

a.  Fan- Compressor  Maps 

The  fan  and  compressor  maps  are  very  similar  and  are  plots  of  corrected  airflow  versus 
pressure  ratio  with  constant  corrected  speed  lines  and  constant  efficiency  islands  (see 
Figure  2).  Entry  to  the  map  is  through  the  corrected  speed  and  Z,  where  Z  is  s  pressure- 
ratio  ratio,  and  is  defined  at  a  Constant  corrected  speed  as  shown  in  Figure  X.  K  is  advan¬ 
tageous  to  use  Z  instead  of  pressure  ratio  because  Z  is  restrained  between  the  limits  of 
0  and  1,  whereas  the  limits  on  the  pressure  ratio  vary  depending  qpon  map  location  and  the 
particular  map.  Also,  an  indication  that  the  fan  or  compressor  is  approaching  surge  is  given 
as  Z  approaches  L. 

b.  Combustor  Map 

The  combustor  map  is  a  plot  of  temperature  rise  across  fibs  combustor  versus  efficiency 
far  constant  input  pressure  (see  Figure  3).  Entry  to  the  map  is  through  temperature  rise  and 
input  pressure,  with  efficiency  being  output. 

c.  Turbine  Map 

The  turbine  map  is  a  plot  of  turbine  corrected  speed  versus  work  function  with  constant 
turbine  flow  fraction  lines  and  constant  efficiency  islands  (see  Figure  4).  The  work  function  and 
flow  function  are  defined  as 
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Figure  2.  Example  of  Fan-Compressor  Map 
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Figure  3.  Example  of  Combustor  Map 


Figure  4.  Example  of  Turbine  Map 
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and 


Entry  to  the  map  Is  through  corrected  speed  and  turbine  flow  function,  with  the  work  function 
and  efficiency  being  output. 

The  work  function  could  have  been  used  as  an  entry  in  place  of  one  of  the  present  entries,  bat, 
because  of  the  shape  of  the  curses,  this  could  lead  to  double  entry  points  far  one  work  ftnctfoo. 
However,  if  the  turbine  maps  were  plotted  in  a  different  format,  this  could  be  an  acceptable 
method. 

2.  DESIGN  POINT 

Once  the  component  maps  have  been  reduced  to  Block  Data  form  and  placed  in  Hie  program, 
it  is  necessary  to  ran  a  design  point.  The  design  point  is  run  at  tboee  conditions  under  which 
the  real  engine  is  designed  or  sized,  usually  sea  level  static.  Design  parameters  necessary  to 
simulate  the  real  engine  (far  example,  airflow,  bypass  ratio,  main  burner  temperature,  various 
pressure  losses,  pressure  ratios,  etc.)  are  input  and  a  complete  thermodynamic  cycle  cal¬ 
culation  is  performed.  For  more  details  on  the  cycle  calculation  see  Section  IV  4,  “Off-Design 
Points.1*  Scale  factors  for  the  component  maps  are  calculated  to  insure  that  the  input  design 
parameters  are  met.  If  the  design  parameters  have  been  correctly  input,  the  design  point  will 
be  completed  after  one  pass  through  the  engine  calculations  (that  is,  no  balancing  will  occur) 
because  Hie  maps  are  shifted  to  reduce  the  errors  to  zero. 

Other  parameters  calculated  and  outptk  at  the  design  point  include  certain  temperatures  and 
airflows,  gas  mixing  areas,  and  nozzle  throat  and  exit  areas. 


3.  SCALING  FACTORS 

Scaling  or  correction  factors  are  calculated  at  the  design  point  using  the  fallowing  equation: 

P  (correction  factor)  =  P  (design)  /  P(map) 

where  P  represents  a  general  parameter.  One  exception  to  this  equation  is  the  equation  for 
calculating  fan  and  compressor  pressure  correction  factors: 

P R (correction  factor)  =  |>R(design)-l  ]/[?**  C«nep)-l] 

where  PR  represents  a  general  pressure  ratio. 

Theoretically,  if  the  component  maps  and  the  input  design  parameters  are  exact  representa¬ 
tions  of  a  particular  engine,  the  correction  factors  will  equal  1.  However,  this  will  not  be  true 
due  to  map  interpolations,  certain  assumptions  such  as  ideal  and  isentropic  How,  and  tolerances 
in  Hie  thermodynamic  calculations.  The  correction  factors  should  be  within  1%  of  JL  Naturally, 
if  unmatched  component  maps  are  used,  the  correction  factors  can  differ  significantly  from  1. 

4.  OFF-DESEN  POINTS 

The  following  discussion  pertains  particularly  to  off-design  points,  alihotgh  the  input  and  the 
general  cycle  calculations  are  the  same  for  the  design  point.  Throughout  the  following  dis¬ 
cussion,  it  should  be  remembered  Hat  scaling  or  correction  factors  (multipliers)  are  to 

all  performance  maps  (Block  Data  parameters). 


AFAPL-TR- 67-125 
Part  I 

For  more  detailed  information  on  the  thermodynamic  equations  used  throughout  the  cycle 
calculation*,  see  References  1,  2,  and  3  and  Part  Hof  this  report.  A  schematic  diagram  of  the 
engine  components  and  station  designations  is  shown  in  Figure  5. 

a.  Input 

The  program  uses  eooutroUed output;  that  is,  the  variables  desired  as  output  can  be  selected 
at  the  start  of  run.  This  selection  is  obtained  by  placing  the  names  of  the  variables  in  the 
first  section  cl  Input  cards.  Controls,  scaling  or  correction  factors,  and  operating  conditions 
make  ip  the  rest  of  the  input. 

The  control  izputs  are  used  to  determine  die  type  of  engine:  mixed  flow  or  separate  flow, 
afterburning,  duct  burning,  and  convergent  or  convergent-divergent  nozzle.  The  controls  are 
also  need  to  fix  the  mode  of  operation:  constant  PCNC,  constant  T4,  or  constant  WFB.  Other 
oontrols  determine  Inlet  conditions,  title  printout,  and  cycle  looping  printouts.  The  correction 
factors  can  be  input  directly,  or  the  design  point  can  be  run  first  and  the  calculated  factors  will 
be  left  In  common.  The  operating  conditions  include  the  flight  Mach  number.,  altitude,  power 
sett  fag  (either  PCNC,  T4,  or  WFB),  duct  burner  and  afterburner  temperatures  or  fuel  Cows, 
bleed,  and  horsepower  extracted. 

b.  feiitial  Values 

The  program  uses  four  primary  independent  variables:  ZF,  PCNF.  ZC,  and  PCNC  (T4  may 
be  substituted  for  PCNC,  depending  tgton  the  mode  of  operation).  Two  secondary  independent 
variables  (TFFHP  and  TFFLP)  are  also  used  to  insure  correct  entry  into  the  turbine  maps. 
Initial  values  for  these  six  variables  must  be  obtained  to  start  the  program  at  each  point.  A 
subroutine  supplies  these  variables  as  a  function  of  T2,  T21,  and  some  of  the  variables  them¬ 
selves.  ft  is  important  to  note  that  the  closer  the  initial  values  are  to  the  final  values  at  a 
balanced  point,  the  faster  the  program  will  run.  Therefore,  after  a  particular  engine  con¬ 
figuration  has  been  rtn  a  few  times,  it  is  usually  advisable  to  change  the  general  initial  value 
equations  to  suit  the  engine,  using  the  knowledge  gained  from  past  runs  to  estimate  more 
closely  toe  final  values  of  toe  variables. 

c.  inlet 

The  thermodynamic  properties  of  toe  atmosphere  are  found  from  a  1952  ARDC  Atmosphere 
Tables  subroutine.  Using  conservation  of  energy  and  issrtroplc  flow,  the  conditions  at  the  face 
of  Oe  fan  can  be  found.  A  ram  recovery  (total  pressure  recovery)  can  be  input  or,  if  not  input,  a 
ram  recovery  defined  by  MI1-E-5008B  Specifications  will  be  used.  S  desired,  a  T2-P2  direct 
input  mode  is  available,  as  are  provisions  for  nonstandard  day  conditions. 

d.  Fan  and  Compressor 

Block  Data  is  used  to  determine  the  performance  characteristics  of  toe  fan  and  compressor. 
When  Z  and  PCN  are  known,  toe  pressure  ratio,  corrected  airflow,  and  efficiency  can  be  found 
by  using  a  general  Block  Data  interpolation  routine  named  SEARCH.  With  the  pressure  ratio 
known  and  when  the  assumption  of  isentropic  compression  and  the  efficiency  are  used,  toe 
thermodynamic  conditions  at  the  exit  of  both  the  fan  and  the  compressor  can  be  calculated. 
Bleed  for  consumer  use,  leakage,  or  cooling  is  accounted  for.  Actual  airflow  leaving  toe  fan  and 
toe  compressor  is  calculated  from  the  corrected  airflow,  temperature,  pressure,  and  bleed. 

e.  Combustor 

The  pressure  drop  in  toe  combustor  is  a  function  of  a  design  pressure  drop  and  toe  ratio 
of  corrected  airflow  to  toe  design  corrected  airflow.  Combustor  efficiency  is  obtained  from 
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Block  Data  using  SEARCH.  The  fuel  used  is  assumed  to  be  JP-4  (at  59*F),  and.  with  the 
assumption  of  adiabatic  and  constant  pressure  combustion,  a  fuel  beating  value  equation  as 
a  function  of  T4  has  been  derived.  Thus  the  fuel/air  ratio,  fuel  flow,  and  thermodynamic  con¬ 
ditions  at  the  combustor  exit  can  be  calculated.  If  WFB  is  known  instead  of  T4,  a  small 

Jtaratibb  is  neoessary. 

f.  Turbines 

Both  turbine  subroutines  use  similar  logic  and  obtain  their  performance  characteristics 
from  Block  Data  using  subroutine  SEARCH.  All  three  turbine  parameters  (CN,  TFF,  DHTC) 
can  be  calculated  before  entering  the  turbine  map,  but  only  two  are  needed.  Therefore,  the 
third  parameter  obtained  from  the  map  is  compared  with  the  calculated  third  parameter,  and  a 
balancing  error  is  generated  if  they  are  not  equal,  hi  this  program,  CN  and  TFF  are  used  for 
map  entries,  and  DHTC  is  used  to  generate  the  error.  In  addition,  the  efficiency  is  also  obtained 
through  SEARCH. 

la  addition,  another  error  will  be  generated  if  TFF  is  not  within  map  limits.  The  error  will 
be  die  difference  between  TFF  and  the  nearest  map  limit.  This  error  becomes  particularly 
important  when  die  estimated  initial  values  of  the  independent  variables  are  far  from  the 
correct  values,  and  the  point  is  extremjly  unbalanced.  When  either  TFF  or  CN  is  not  within 
map  limits,  they  are  set  to  the  nearest  map  limit,  and  one  of  Hie  independent  variables  is 
changed  in  an  attempt  to  rectify  Hie  situation.  The  operating  point  must  appear  on  all  maps 
before  a  complete  cycle  calculation  can  be  accomplished. 

Horsepower  extraction  is  accounted  for  in  calculating  DHTC  of  Hie  high  pressure  turbine. 
When  Hie  efficiency  is  used  and  Hie  turbine  process  is  assumed  isentropic.  the  thermodynamic 
propertleu  at  both  turbine  exits  can  be  calculated.  Any  bleed  airflow  for  -iooling  Hie  turbines  is 
treated  as  if  it  entered  the  main  stream  behind  the  turbine,  and  the  thermodynamic  properties 
at  the  turbine  exits  are  recalculated  *o  account  for  this. 

g.  Duct 

The  duct  airflow  and  bypass  ratio  are  calculated  from  the  fan  and  compressor  airflows.  The 
pressure  drop  in  Hie  duct  is  treated  e®  in  the  main  combustor.  For  duct-burning,  the  same 
fuel  heating  value  equation  that  was  used  in  the  combustor  is  again  used,  but  the  efficiency  must 
be  iqput.  As  in  Hie  combustor,  either  the  temperature  (T24)  or  the  fuel  flow  (WFDj  may  be  input. 

If  a  separate  flow  engine  is  being  simulated,  the  duct  nozzle  calculations  are  done  in  this 
routine,  although  they  are  accomplished  in  the  same  manner  as  for  the  main  nozzle. 

h.  Mixer 

The  gas  mixing  areas  (duct  exit  and  turbine  discharge  for  a  mixed  flow  engine  or  just  the 
turbine  discharge  area  for  a  separate  flow  engine  are  calculated  at  Hie  design  point  using 
either  an  input  static  pressure  or  Mach  number.  At  an  off-dcsign  point  the  areas  are  used  to 
calculate  static  pressures  and  Mach  cumbers. 

Far  a  separate  flow  engine.  Hie  thermodynamic  conditions  entering  the  afterburner  are  cow 
known,  since  they  are  identical  to  turbine  discharge  conditions. 

Far  a  mixed  flow  engine,  a  set  cl  derived  equations  based  cm  one-dimensional  fluid  flow 
theory  and  conservation  of  mess,  energy,  and  momentum  is  used  to  determine  the  thermo¬ 
dynamic  conditions  after  complete  mixing  of  the  two  gas  streams  (Reference  4J,  These 
equations  do  not  require  that  the  static  pressures  of  the  two  entering  streams  be  equal.  How¬ 
ever,  for  a  correct  engine  balance,  the  two  static  pressures  must  be  equal,  and  a  bciaccing 
error  is  generated  if  they  are  not  equal. 
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i.  Afterburner 

The  dry  loss  (cold  loss)  pressure  drop  in  the  afterburner  is  a  function  of  £  design  presauro 
drop  and  the  ratio  of  corrected  gas  flow  to  the  design  corrected  gas  flow. 

For  afterburning,  the  same  equation  for  the  fuel  heating  value  that  was  used  in  the  combustor 
is  again  used,  but  the  efficiency  must  be  input.  As  in  the  combtstor,  either  the  temperature 
(T7)  or  the  fuel  flow  (WFA)  may  be  input.  A  momentum  loss  (hot  loss)  pressure  drop  is  also 
calculated. 

j.  Nozzle 

The  main  nozz*  2  program  uses  fixed  effective  areas  (except  when  afterburning)  calculated 
at  the  design  point.  Either  a  convergent  or  a  convergent-divergent  subroutine  may  be  used 
depending  upon  the  input  controls.  If  afterburning  has  been  selected,  the  nozzle  areas  are 
allowed  to  float  to  obtain  optimum  performance:  however  the  areas  are  returned  to  their 
original  design  values  after  the  afterburning  point  is  completed.  The  duct  nozzle  behaves 
identically  with  the  main  nozzle,  including  floating  areas  if  duct-burning  has  been  selected. 

Because  all  thermodynamic  properties  of  the  gas  stream  are  known,  as  well  as  the  amount 
of  flow,  nozzle  areas,  and  ambient  pressure,  there  is  a  redundant  parameter.  For  this  program, 
the  total  pressure  of  the  gas  stream  was  chosen  as  the  redundant  parameter.  The  nozzle  cal¬ 
culations  (Reference  5)  are  made  without  using  the  total  pressure,  and  a  required  total  pressure 
compatible  with  all  other  known  parameters  is  calculated.  This  required  pressure  is  compared 
with  the  actual  pressure,  and  a  balancing  error  is  generated  if  they  are  not  equal. 

k.  Performance  and  Output 

At  this  point,  six  errors  have  been  generated  after  one  pass  through  the  engine.  Several 
more  passes  must  be  completed  under  control  of  the  error  matrix  and  engine  balancing  sub¬ 
routines.  See  Section  V  for  a  detailed  description  of  the  balancing  technique  used.  Eventually, 
however,  the  errors  wail  be  reduced  to  zero,  and  engine  performance  will  be  calculated  using 
standard  equations.  Gross  thrust  is  obtained  by  summing  the  momentum  term  xa  nozzle 
velocity  coefficient  may  be  input)  and  pressure-area  term,  and  net  thrust  is  in  turn  found  by 
subtracting  a  ram  drag  (airflow  momentum  loss  at  inlet)  term  from  the  gross  thrust.  Specific 
fuel  consumption  is  total  fuel  flow  divided  by  net  thrust. 

As  previously  mentioned,  a  controlled  output  is  used,  whereby  only  selected  variables  are 
printed.  Each  variable  is  labeled  with  its  name,  and  provisions  have  been  made  for  changing 
the  name  of  a  variable,  hi  addition,  the  values  of  all  variables  in  common  are  printed  in  a  close 
format  so  that  variables  other  than  those  selected  for  a  specific  run  are  available  later  on. 

5.  QUADRATIC  INTERPOLATION  ROUTINE 

Throughout  foe  program  there  are  many  small  loops  (far  example,  thermodynamic  iterations 
and  table  look-ip)  which  require  convergence.  Trial- and-error  methods  and  linear  interpola¬ 
tions  can  be  lime-consuming,  especially  when  a  tight  tolerance  is  necessary;  therefore  a 
general  interpolation  routine  called  AFQUIR  (Air  Force  Quadratic  Interpolation  Routine)  was 
developed. 

This  routine  requires  a  dummy  array  dimensioned  for  nine  locations.  Also  input  into  the 
routine  through  the  calling  argument  axe  the  independent  and  the  dependent  variables,  the 
answer  or  value  which  the  dependent  variable  is  to  converge  qpon,  the  number  of  tries  at  ooo? 
vergence,  foe  tolerance,  and  a  variable  called  DIR. 
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The  DIR  is  either  set  or  calculated  in  the  calling  program  and  is  an  initial  guess  at  the 
direction  and  percentage  change  to  apply  to  the  first  -value  of  the  independent  variable.  If  not 
mortal  is  known  about  the  variables  to  calculate  a  DIR,  an  arbitrary  value  may  be  set.  This 
should  not  affect  the  final  result,  but  may  increase  the  cumber  of  tries  at  convergence. 

TbeDIR  thus  establishes  the  second  value  of  the  independent  variable.  This  value  is  used 
is  the  calling  program  to  determine  a  corresponding  second  value  of  the  dependent  variable 
sad  ATQUIR  is  called  a  second  time  with  two  sets  of  values.  A  linear  interpolation  is  made 
which  results  in  a  third  value  nf  the  independent  variable.  AFQUIR  is  then  called  a  third  time 
with  the  third  values  of  independent  and  dependent  variables  and  a  quadratic  interpolation  is 
made.  The  values  of  time  three  sets  of  variables  have  been  stored  in  the  dummy  array,  and 
from  here  on,  quadratic  interpolations  are  made  using  the  three  sets  which  give  values  closest 
to  *he  answer.  Values  farthest  from  the  answer  are  lost. 

Various  safeguards  are  built  into  AFQUIR  to  return  the  interpolation  method  to  DIR  or 
linear  if  the  roots  of  the  quadratic  become  complex,  if  the  quadratic  does  not  intercept  the 
anewor,  if  the  value  of  the  independent  variable  differs  radically  from  previous  values,  or  if 
two  sets  of  independent  and  dependent  variables  are  vy^tieai. 

Also,  it  is  possible  to  preload  the  dummy  array  and  to  start  directly  at  the  linear  or 
quadratic  interpolations  if  desired. 

hi  summary,  AFQUIR  is  aoompletely  flexible  routine  which  performs  quadratic  interpolation 
for  quick  convergence  ol  general  functions. 
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SECTION  V 

BALANCING  TECHNIQUE 

The  balancing  technique  is  based  upon  finding  a  solution  for  a  set  of  partial  differential 
equations.  For  this  program,  the  set  is  composed  of  six  equations;  however,  using  a  set  of 
only  three  equations  will  simplify  the  follcrwing discussion.  This  corresponds  to  a  basic  turbo¬ 
jet  engine  simulation.  It  is  relatively  easy  to  expand  the  set  of  three  equations  to  one  of  six, 
as  required  in  SMOTE,  or  even  further.  For  example,  a  triple-spool  turbofac  would  require 
nine  equations. 

As  discussed  previously,  six  independent  variables  were  selected  (ZF,  PCNF,  ZC9  PCNC  or 
T4,  TFFHP,  and  TFFLP).  Once  these  variables  have  been  given  initial  values,  it  is  possible 
to  proceed  through  an  entire  engine  cycle  calculation.  Six  errors  are  generated  as  shown  in 
Section  IV.  These  initial  values  of  the  six  variables  and  six  errors  are  referred  to  as  base 
values. 

In  the  following  equations,  V  refers  to  a  variable  and  E  to  an  error.  The  basic  set  of  differ¬ 
ential  equations  based  on  E  =  f  (V)  is  (Reference  6) 
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where  the  single  subscripts  correspond  to  three  variables  and  three  errors  and  where  the 
double  subscripts  indicate  the  change  in  a  particular  error  (first  subscript)  due  to  a  change  in 
a  particular  variable  (second  subscript). 

Assuming  small  changes  results  in  the  following  approximations  (where  B  refers  to  a  base 
value): 
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With  these  approximations  a- jd  the  fact  that  E  should  be  zero  when  the  engine  is  balanced, 
the  set  of  partial  differential  equations  reduces  to 


E,  -  EB, 


e*-eb2 


e3-eb. 


ae„ 

dV, 

A  . 

AE,* 

<r^  +  - 

AEa 

dV,  = 

-EB, 

Av, 

T 

AV* 

AV3 

AE*, 

Av, 

dV, 

+ 

AE** 

AV* 

dV*-f 

AE*3 

AV, 

dV,= 

-EB* 

ae3, 

dV, 

i  - 

AE3* 

dV*  +- 

ae33 

dV 

"EB, 

AV, 

T 

Av2 

AVs 

13 


V 


AFAPL-TR-67-125 
Part  I 

Three  more  passes  {sis  for  SMOTE)  ”»*e  now  made  through  the  engine  cycle  calculations, 
and  one  variable  is  changed  by  a  small  amount  ( AV )  for  each  pass.  'P'e  charge  in  each  error 
due  to  the  small  change  in  the  variables  (Ae/AV)  can  then  be  calc  o’  _edL 

The  above  set  of  differential  equations  can  now  be  solved  for  dVy,  dVg,  and  dVg,  and,  in 
general,  the  new  value  of  each  independent  variable  would  be  given  by 

V  =  VB  +  dV 

If  the  engine  cycle  calculations  were  linear  functions,  the  engine  would  balance  {errors  equal 
zero)  wife  these  new  values  of  fee  variables.  However,  this  is  not  fee  case,  and  it  is  usually 
necessary  to  repeat  fee  above  process  (where  fee  new  values  become  fee  base  values)  several 
times  before  a  balance  is  obtained. 

A  subroutine  to  determine  the  solution  of  a  matrix  is  used  to  solve  fee  set  of  differential 
equations.  After  each  pass  through  fee  engine,  a  matrix  array  is  loaded  wife  fee  appropriate 
values;  after  seven  passes  (base  value  plus  six  independent  variables),  the  matrix  subroutine 
is  called  to  solve  the  matrix. 

£  was  found  feat  the  s‘dV*s”  obtained  from  the  solution  of  fee  differential  equations  were  in 
many  cases  too  large,  feus  causing  the  variables  to  exceed  their  limits,  ana  to  make  it  practi¬ 
cally  Impossible  to  balance  the  cycle.  The  “dV*s”  are  therefore  multiplied  by  a  suppression 
factor  (presently  0.6)  which  limits  fee  swing  of  fee  variables.  la  addition,  if  a  suppressed 
**dV”  is  still  greater  than  5%  of  the  value  of  fee  variable  itself,  it  is  reduced  to  fee  5%  value. 
Although  this  procedure  may  tend  to  increase  the  number  of  passes  before  balancing  in  some 
cases,  it  also  balances  points  which  previously  would  not  balance.  These  points  are  most 
generally  far  from  fee  design  point,  where  oscillations  of  fee  dependent  variables  tend  to  build 
vp. 
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APPENDIX 


SAMPLE  RESULTS 


The  following  computer  printouts  are  examples  of  typical  output  from  SMOTE.  The  first 
point  is  tbs  design  point  and  includes  a  page  of  correction  (or  scaling)  factors  and  a  page  of 
values  of  variables  in  common.  The  other  points  represent  conditions  throughout  a  flight 
envelope  and  consist  of  a  primary  page  of  output  for  each  point.  Not  Included  fur  these  points 
is  a  common  dump,  which  normally  follows  each  primary  output  page  and  is  very  similar  to 
the  common  dump  following  the  design  point  correction  factors. 

The  engine  cycle  chosen  was  a  mixed  flow  turbofan  (bypass  ratio  of  1.4)  with  a  convergent 
nozzle,  a  total  airflow  of  180  pounds  a  second,  and  a  turbine  inlet  temperature  (T4)  of  2400*11. 
The  points  were  run  in  a  fixed  T4  mode;  that  is,  PCNC  is  an  independent  variable.  No*®  that 
the  coszle  area  is  recalculated  at  each  afterburning  point  for  optimum  expansion  and  that  no 
balancing  occurs  at  these  points. 
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i-  This  report  describes  a  digital  computer  program  titled  SMGTE  (Simulation  of  Turbofan 
Engine).  SMOTE  is  a  computer  program  for  balancing-cycle  tuxbofsn  engines  capable  of 
running  both  design  and  off-design  points.  Component  performance  maps  are  reduced  to 
Block  Data  (tabular  form)  to  provide  a  base  for  calculating  component  performance.  The 
design  point  is  ran  first  and  map  correction  factors  are  calculated  to  scale  the  components 
to  the  desired  performance.  These  correct! on  factors  are  then  applied  to  the  component 
performance  maps  at  off-design  points.  Initially,  when  tbs  program  Is  running  at  an  off- 
design  point,  the  cycle  is  cot  balanced-  and  errors  (for  example,  work  required  by  the 
compressor  miens  work  supplied  by  the  turbine)  are  generated.  Small  changes  in  engine 
independent  variables  (for  example,  compressor  speed)  then  produce  small  changes  in  the 
errors,  and  these  differential  changes  dre  loaded  into  &  matrix.  The  matrix  is  then  sorted 
for  the  set  of  Independent  variables  which  results  in  zero  errors,  thus  balancing  the  cycle. 
Actually,  this  process  may  be  repeated  several  times  before  it  reaches  a  balanced  point 
because  there  is  a  non  linear  relationship  between  lie  independent  variables  and  the  errors. 
Sample  results  arc  included  in  this  report. 

(Distribution  of  this  abstract  is  unlimited.) 
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